Angular Dependence of the Nonlinear Transverse Magnetic Moment of YBa 2 Cu 3 06.95 

in the Meissner State. 
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The angular dependence of the nonlinear transverse mag- 
netic moment of untwinned high-quality single crystals of 
YBa2Cu3C>6.95 has been studied at a temperature of 2.5K 
using a low frequency AC technique. The absence of any 
signature at angular period 2ir/4 is analyzed in light of the 
numerical predictions of such a signal for a pure d x 2_ y 2 order 
parameter with line nodes. Implications of this null result 
for the existence of a non-zero gap at all angles on the Fermi 
surface are discussed. 



Measurements of the low energy quasiparticle exci- 
tation spectrum of high temperature superconductors 
(HTS) are crucial in elucidating the nature of the pair- 
ing state and gap function in these materials. For a 
d x 2_ y 2 order parameter, the gap function varies as A(tp) 
= \Aocos(2ip) |, riding on a cylindrical Fermi surface in 
momentum space, with zeros or line nodes at (p = (tt/4 + 
rm/2). Here, A is about 25meV |j] for YBa 2 Cu 3 6 .95 
(YBCO), and <p is the angle measured from the crystallo- 
graphic a or 6 directions. The density of states for quasi- 
particle excitations near these nodes on the Fermi surface 
increases linearly with energy. In a two fluid model, this 
results in the superconducting condensate density n s (T) 
decreasing linearly with increasing temperature from its 
zero temperature value n s (0) ||. In conventional s-wave 
superconductors, there is a non-zero gap everywhere on 
the Fermi surface for quasiparticle excitations, and the 
depletion of the condensate is exponentially small at the 
lowest temperatures 

The temperature dependence of the in-plane London 
penetration depth || X(T) in these materials has been 
widely cited as evidence in support of the presence of line 
nodes in the gap function, though at the lowest temper- 
atures there is some deviation from the expected linear 
behavior, even in the best samples Measurements 
of the magnetic field dependence of the low temperature 
specific heat have also been interpreted as evidence for 
the presence of nodes, but the analysis is subject to the 
limitations of a many-parameter fit ||. Data from later 
experiments, taken over a wider range of temperature 
and field do not yield the same conclusions upon anal- 
ysis Q. Recent studies on the scaling of specific heat 
data has also been cited as evidence in support of nodes 
, but the scaling works poorly, and only over a limited 



temperature range. Amongst experiments that provide 
angular information, angle resolved photemission exper- 
iments (ARPES) can be interpreted as showing a mini- 
mum in the energy gap in the (110) direction in YBCO 
and BSCCO (l]]§], but have a resolution of only a few 
meV (lmeV = 11.62K) and thus cannot resolve excita- 
tions or an underlying s-wave gap at the lowest energies. 
Inelastic neutron scattering shows evidence in support of 
minima ||, but there too, the resolution is lmeV at best. 
Thus, one can safely say that the overall picture is quite 
ambiguous regarding the possibility of an underlying gap 
below lmeV. 

In this Letter, we describe an experimental technique 
for distinguishing nodes from deep minima (quasinodes) 
in the order parameter, that also provides information 
about the angular position of these nodes or quasin- 
odes, by probing the existence of low lying excitations 
in response to an applied magnetic field in the Meissner 
regime. Here, the kinetic energy of the superflow of the 
screening currents provides the energy for quasiparticle 
excitations. Our null result for this probe rules out the 
existence of nodes, and allowing for quasinodes, sets a 
lower bound on the size of the underlying gap. 

For a type II superconductor in a magnetic field in 
the Meissner regime, screening supercurrents flow in a 
volume near the surface given approximately by the pen- 
etration depth A. For the condensate participating in 
these currents, the quasiparticle excitation spectrum is 
modified by a semiclassical 'Doppler shift' to E(k) — 
y/ (A| + e|) + v s ■ vp 10 1, where v s is the superfluid 
flow field and vp is the Fermi velocity. For a gap function 
with nodes, this leads to quasiparticle excitations even at 
zero temperature @ . Due to the linear density of states 
near a line node, the depletion of the condensate due to 
quasiparticles is proportional to (v s -vp) 2 . These quasi- 
particles create a 'backflow' which is then responsible for 
a nonlinear contribution to the magnetization that goes 
as FP\ 2 , where H is the applied magnetic field, and v s 
~ HX. In the case of a d x 2_ y 2 order parameter, the gap 
function has a four-fold angular symmetry, and this gives 
rise to an intrinsic nonlinear transverse magnetization, 
superimposed on the nominal diamagnetic response |L2| . 
Since this effect is felt only by a fraction of the conden- 
sate within a thickness of A of the surface, the nonlinear 
magnetic moment is proportional to the surface area of 
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the sample, and not the full sample volume. Also, since 
A >>£, where £ is the in-plane coherence length, this is 
a bulk effect. 

We are interested in the nature of the gap function 
for supercurrents flowing within the Cu-0 planes of the 
superconductor, which consists of layers of these planes. 
The YBCO crystals are flat with the c-axis oriented per- 
pendicular to the crystal plane which is the a-b plane. 
The magnetic field is applied parallel to the crystal plane 
and the transverse magnetic moment is measured in the 
crystal plane perpendicular to the applied field. As the 
crystal is rotated in the applied magnetic field about the 
c-axis, the screening currents flow in different directions 
relative to the in-plane gap function, which is pinned to 
the crystal lattice. This leads to an angular modula- 
tion of the nonlinear transverse magnetic moment that 
provides angle-resolved information about the low lying 
excitations in YBCO. 

Numerical calculations of the nonlinear transverse 
magnetic moment (m?) of a sample with a finite disk 
shaped geometry have been carried out for pure d-wave 
and mixed order parameters. For pure d-wave, the cal- 
culation predicts the amplitude of the expected four-fold 
modulation of tut with angular period n/ 2 p"2| ], a conse- 
quence of the symmetry of the gap and angular position 
of the nodes. The presence of a small s-wave component 
as in d+s changes the angular position of the line nodes, 
and introduces a component in the angular modulation 
of m,T with period n p3[ . However, the angular period 
7r / '2 component is not adversely affected by small ad- 
ditions of s-wave. Anisotropy of the Fermi surface and 
anisotropy in A a and A& have also been considered, and 
have very similar consequences. For experimentally de- 
termined values of such anisotropy jyj, there is a small 
additional angular period 7r component, but the period 
ir/2 component is essentially unchanged. On the other 
hand, for a d-wave like order parameter, with varying lev- 
els of a non-zero quasinode, as in d x 2_ y 2 + is and d x i_ y i + 
id xy symmetries, the amplitude of the period n/ 2 com- 
ponent is suppressed (see Fig. 2 inset). In addition, to 
prevent the nonlinear Meissner effect (NLME) from be- 
ing thermally washed out, the temperature has to be such 
that where H a = ■^tfz, and <j> is the flux quan- 

tum. For \ ab = 1400A and £ a b = 20 A, which are typical 
values, H « 8000 Oe. For measurements in a field of 
amplitude 300 Oe, this yields jy- = 0.0375, and requires 
that T < 10K. 

To measure mr, an AC magnetic field is applied in the 
a-b plane of a mm size single crystal of YBCO. The field is 
modulated at 12Hz, and the transverse magnetic moment 
at 36Hz is detected by the transverse coil of a Quantum 
Design SQUID susceptometer. The signal is expected at 
the third harmonic because hit ~ sgn(H)lP. The sample 
is maintained at 2.5K with continuous cooling through- 
out the measurements. The superconducting magnet is 
operated with its persistent-current switch open, and the 
AC magnetic field is generated by driving a very low dis- 



tortion (<-110dB in higher harmonics) and low noise AC 
current through the magnet coil, with a maximum ampli- 
tude of 300 Oe. The sample is placed in the center of the 
magnet coil and aligned optimally with the transverse 
flux detection coils. The analog output from the trans- 
verse SQUID is fed directly to a phase sensitive detection 
system. The large 12Hz background from the fringing 
field of the magnet and geometric demagnetization fields 
from the sample are rejected by a high-pass filter section. 
The signal at 36Hz is then detected with a phase sensi- 
tive detector locked to Sf of the current generator. The 
sample is rotated in situ between measurements in steps 
of 6° using a custom built sample holder that allows for 
in situ detection of the angular position of the sample 
at low temperatures. This system has been described in 
detail elsewhere | jl5"| ]. It allows for precision in angular 
positioning of about 0.1° between steps and accuracy of 
better than 1° in an entire rotation. The sample holder 
and associated angular detection systems are designed to 
provide minimal magnetic background. 

The detection setup has been calibrated and tested by 
'simulating' a real magnetic moment in an environment 
identical to that of the actual experiment. This is done 
by running an AC current through a copper coil mounted 
on the sample stage at a frequency 5/, referred to the fre- 
quency /=12 Hz of the oscillating magnetic field, which 
has an amplitude of 300 Oe. The AC moment of 1.5 x 
10 -8 emu rms amplitude at 5/ is detected in the presence 
of the large background signal at frequency / due to the 
fringing field of the magnet. This calibration was done 
at 2.5K, with just the bare sample holder. The measured 
amplitude is shown in Fig.l as a function of the angular 
orientation of the coil. The angular Fourier transform of 
this data (inset) shows the amplitude at angular period 
2ir. The 'noise floor' in the Fourier transform of the cal- 
ibration signal at higher angular frequencies corresponds 
to an amplitude less than 5 x 10~ 10 emu. 

The angular dependence of the NLME has been mea- 
sured in a number of untwinned single crystals of YBCO, 
including a very high quality crystal (UWHC), grown in a 
YSZ crucible by the Argonne group. This crystal, which 
is 99.95% pure, has been cut into a disk shape with diam- 
eter 1.5mm, and is 67 fim thick. The high quality of the 
crystal is borne out by a rocking curve with full width 
at half maximum of 0.08° for the 006 peak of a high res- 
olution x-ray scan, and there were no twins observable 
in repeated area scans over different parts of the crystal. 
The superconducting transition of this crystal has an on- 
set of 93K and a width of less than 2.5K as measured 
magnetically with a 10 Oe field applied in the a-b plane. 
The field of first flux entry was measured to be about 300 
Oe in the a-b plane at 2.5K. 

Measurements of tut are shown in Fig. 2 as a function 
of the angular orientation of the sample with respect to 
the applied magnetic field, with the 'a' direction being 
initially oriented along the field. The sine Fourier trans- 
form of the data is shown in the inset. For a pure d- 
wave order parameter, we are interested in the angular 
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period 2it/4 component of nvr, the predicted value for 
which is 1.7 x 10 _9 emu according to the calculations in 
Ref. |l^] . This component is clearly in the noise of the 
data, below 5 x 10~ 10 emu. This 'noise' is in part due to 
residual trapped flux, and also due to the noise floor of 
the measurement apparatus, as has been verified in re- 
peated measurements. This imposes an upper bound on 
the size of the nonlinear Meissner effect, which if there, is 
less than 30% of the predicted value. Identical measure- 
ments have also been carried out on a very high quality 
rectangular crystal (UBCyca) grown by the UBC group 
in a BaSZ crucible It has the same surface area as 
our disk shaped crystal (UWHC) and the results were 
essentially the same. 

This null result has to be examined in the light of mea- 
surements of the penetration depth. Scrutiny of data 
on the temperature dependence of the penetration depth 
X(T) indicates that there is always some curvature away 
from the pure d-wave result at the lowest temperatures. 
This is found in all published data on the cleanest bulk 
single crystals , and has been attributed to the pres- 
ence of unitary scattering centers, that do not affect 
T c |l^, |l7). However, even with the concentration of 
unitary scatterers needed to produce the required cur- 
vature, the nonlinear transverse magnetic moment effect 
may only be reduced to about 90% of the full value jll) , 
and our data indicates a much stronger suppression. This 
same data for \(T) for UBCyca can be fit at the low- 
est temperatures by a model which contains quasinodes, 
i.e., where the is 2.5%, A min being the residual gap 
in the nodal direction. Such a fit is significantly bet- 
ter at the lowest temperatures than the fit for a pure 
d-wave order parameter. The suppression of the nonlin- 
ear transverse magnetic moment in all our measurements 
is consistent with A £ in between 2 - 3%. In the event that 
the X(T) and photoemision results can be explained by 
something other than a minimum in the gap, our experi- 
ment is also consistent with a d+s order parameter where 
A s >A d . 

Recently, other experiments that have attempted 
Jl9| , p0[ to see the NLME in measurements of the field 
dependence of A at low temperatures, but have not ob- 
tained conclusive results. According to theory for a 
pure d-wave order parameter, A should vary linearly with 
H at the lowest temperatures where ^~ < ~§~ ■ These most 
recent measurements were done at a very high level of 
sensitivity, about 100 times better than an earlier effort 
pl[ . In one of these experiments done at UIUC |l9| , a sig- 
nificant field dependence in A was observed, but this was 
attributed to trapped flux as the field dependence was 
closer to \/i?and the temperature dependence was not 
as p redicted by theory. The experimental effort at UBC 
ppf observed a field dependence in A, and also measured 
the effect as a function of the in-plane angular orientation 
of the applied field. However, they report that the field 
dependence measured at ± 45° to the crystal axes were 
not identical as would have been expected from theory, 



and the temperature and field dependence of the effect 
could not be understood within the current picture for 
the NLME 0. 

Earlier attempts at measuring the angular dependence 
of the nonlinear transverse magnetic moment using a 
DC technique p2| were less sensitive by more than two 
orders of magnitude than the present work due to the 
very large linear background in a DC measurement, and 
the analysis overestimated the size of the effect signifi- 
cantly. In our experiments, trapped flux may produce 
a signal, but the angular dependence of this effect has 
the largest amplitude at angular period 2tt, and is eas- 
ily distinguished from the period 2~k/4 modulation for 
which we are searching. The angular modulation aris- 
ing from the geometric demagnetization factor is linear 
in the field below the field of first flux entry, and is mini- 
mized by looking at only the nonlinear component which 
is extracted directly in our work. Thus, artifacts due to 
trapped flux and geometry are both minimized by our 
technique. 

There have been some calculations of the effects of non- 
local electrodynamics |^3|,^4j]. These are motivated by 
the idea that in a BCS like picture, when the gap 
goes to zero, then £ = Hvp/irAk diverges, and one might 
no longer be in the local limit where £ < A . However, 
it turns out that these considerations are relevant only 
when H \[c. In our experiment, H is applied in the a-b 
plane, and the volume within a depth A of the a-b plane 
crystal face that is responsible for the nonlinear Meissner 
effect is not affected by nonlocal effects. Even if we con- 
sider a 'weakly 3D' case |24|], the effect is at a field scale 
of about 20 Oe, far too small to suppress the NLME, 
whose characteristic field in our measurements is 300 Oc. 

In summary, in view of the results of this paper, and 
in light of the evidence for a predominantly d-wave order 
parameter |p5| from other experiments, we suggest that 
the order parameter in YBCO may still be d-wave like, 
but have quasinodes instead of line nodes, with / ^ lin be- 
tween 2 - 3%. Measurements of the penetration depth 
to lower temperatures or photoemisssion experiments at 
higher resolution may confirm this. 
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FIG. 1. Sensitivity of the measurement apparatus to a 
magnetic moment of 1.5 x 10 -8 emu at 3f. The amplitude 
at period 2n is shown in the cosine Fourier transform (inset) 
and is 2.2//V. 



FIG. 2. rriT(3f) as a function of angle with with respect to 
applied magnetic field for UWHC YBCO sample. The ini- 
tial orientation is with the 'a' direction parallel to the applied 
magnetic field. The absolute value of the sine Fourier ampli- 
tudes are plotted in the inset. The expected amplitude at 
angular period 2ir/4 modulation for pure d-wave is shown as 
a horizontal line at 17 x 10 -10 emu. The level predicted for 
Amin/ A — 2.5% is shown by the horizontal line at 5 x 10~ 10 
emu, near the noise floor. 
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